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In this short note we clarify the role of the boundary terms in the calculation of the leading order 
tree-level bispectrum in a fairly general minimally coupled single field inflationary model, where the 
inflaton's Lagrangian is a general function of the scalar field and its first derivatives. This includes 
k-inflation, DBI-inflation and standard kinetic term inflation as particular cases. These boundary 
terms appear when simplifying the third order action by using integrations by parts. We perform 
the calculation in the comoving gauge obtaining explicitly all total time derivative interactions and 
show that a priori they cannot be neglected. The final result for the bispectrum is equal to the 
result present in the literature which was obtained using the field redefinition. 

I. INTRODUCTION 

The inflationary scenario has become the most compelling idea for what really happened in the very early stages 
of our universe. The main reason for this is that inflation is very successful at explaining the problems of the old 
Big Bang model and it also provides us with a mechanism to generate the primordial seed perturbations that later 
evolved into the cosmic microwave background radiation anisotropics and the large-scale structure of galaxies. 

The most simple models of inflation predict a nearly Gaussian and nearly scale invariant primordial perturbation 
in good agreement with recent observations. However, a small amount of non-Gaussianity is still allowed by the data. 
If such a small contribution were detected and if it were of primordial origin, it would have profound implications 
for inflationary models and our understanding of the very early universe. Because of this reason, recently there has 
been a huge effort to try to construct models that predict large (i.e. observable) non-Gaussianity and calculate the 
higher-order correlation functions like the bispectrum and the trispectrum. This search has been productive. Many 
possibilities have been found, for example models with non-canonical kinetic terms (like DBI-inflation, k-inflation, 
ghost-inflation) (l|-[l3| , multiple fleld models of inflation p"3l - l43 | , temporary violations of the slow-roll conditions and 
small departures of the initial vacuum state from the standard Bunch-Davies vacuum [iB - fsH - For recent reviews 
about these mechanisms to produce non-Gaussian perturbations, see Ref. (5^ - [56j . 

In this note, we will be interested in the calculation of the bispectrum for quite general models of k-inflation. Our 
conclusions will apply to all k-inflation, DBI-inflation and standard kinetic term inflation models. Our main goal is 
to clarify the role of the boundary terms in the calculation of the bispectrum. These boundary interactions appear 
when one does many integrations by parts to simplify the third order action. 

For the standard kinetic term inflation all these boundary terms have been for the first time recently calculated 
in [syj . However in all calculations of the bispectrum so far these terms have been neglected and a field redefinition was 
used instead. Maldacena ^58^ states that these boundary interactions are important and he takes them into account 
using field redefinitions. Wc should point out one exception 59], where the authors working in the uniform curvature 
gauge and considering standard kinetic term inflation show explicitly that the bispectrum can be determined by using 
the third order action without the need to redefine the field if the contribution from the boundary terms is included. 

In this work, we will perform an analogous calculation, but this time for a general k-inflation model and working in 
the comoving gauge. In this gauge, one works from the beginning with the variable of interest, the comoving curvature 
perturbation. Naively, the third order action does not seem to be suppressed by the slow-roll parameters contrarily 
to the expected. However it has been shown [1, [s^ that after many integrations by parts this suppression becomes 
evident. This procedure produces many boundary terms which are the main focus of this work. 

This paper is organized as follows. In section |TT1 we will introduce the model, the background spacetime and some 
useful notation. In section Hill we will discuss linear perturbations and present the solution for the mode functions. 
In section IIVI we shall review the usual calculation of the third order action and the bispectrum and obtain the new 
boundary terms which we argue can be used to calculate the bispectrum without the need to do a field redefinition. 
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We also present the calculation of the bispectrum produced by these time boundary interactions. Section |V] is devoted 
to the conclusions. 



II. THE MODEL 



In this note, we will consider the class of models described by the following Lagrangian 

S^IJ d^^V^ [MliR + 2P(X, 0)] , (1) 

where (j) is the inflaton field, Mpi is the reduced Planck mass, R is the Ricci scalar, X = —\g^'^d^(t)d^4) is the 
inflaton's kinetic energy and g^^u is the metric tensor. P{X, (f) denotes the inflaton's Lagrangian and we assume it 
is a well-behaved function of its two variables. Throughout this work, we use a system of units where the Planck 
constant /i, the speed of light c and the reduced Planck mass Mpi are set to unity. This general Lagrangian includes 
as particular cases the DBLinflation model [2, 60] and the k-inflation model 61]. 

We are interested in flat, homogeneous and isotropic Friedmann-Lemaitre- Robertson- Walker background universes 
described by the line element 

ds^ = -dt^ + a^{t)5^jdx^dx' , (2) 

where a{t) is the scale factor. The Friedmann equation and the continuity equation read 

iH^^p, p^-3Hip + P), (3) 

where dot denotes derivative with respect to cosmic time, the Hubble rate is H = d/a, p is the energy density of the 
inflaton and it is given by 

p = 2XPx - P, (4) 

where Px denotes the derivative of P with respect to X. It was shown in [g^] that for this model the speed of 
propagation of scalar perturbations ( "speed of sound" ) is Cg given by 

c^^^= (5) 
^ p,x P.x + 2XP,xx ■ ^ ' 

We define the slow- variation parameters, analogues of the slow-roll parameters, as 

H XP^x e c; 

^ = -^ = ^' '=^- 

We assume that the rate of change of the speed of sound is small (as described by s) but Cg is otherwise free to change 
between zero and one. It is convenient to introduce the following parameters that describe the non-linear dependence 
of the Lagrangian on the kinetic energy 

E = XPx + 2X^Pxx = A = X^Pxx + ^X^P.xxx- (7) 

cj 3 

These parameters are related to the size of the bispectrum. 



III. PERTURBATIONS 



In this section we will consider linear perturbations of the background ([2]). There is a vast amount of works on 
linear perturbations, see for example [64] • It is convenient to use the Arnowitt, Deser and Misner (ADM) metric 
formalism [g^. The ADM line element reads 

ds^ = -N'^df + {dx' + N'dt) {dx^ + N^dt) , (8) 

where N is the lapse function, N'' is the shift vector and hij is the 3D metric. 
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The action ([T]) becomes 

S=\j dtd^xVhN + 2P) + i J dtd^xVhN-^ {EijE'^ - E^) . (9) 



The tensor Eij is defined as 



= - % - V, A^, - V,iV, , (10) 



and it is related to the extrinsic curvature of the spatial slices by Kij = N~^Eij. is the covariant differentiation 
with respect to hij and all contra- variant indices in this section are raised with hij unless stated otherwise. 
The Hamiltonian and momentum constraints are respectively 

^^'>R + 2P - 2n^N-^P^x - {E,jE'^ - E^) = 0, 

Vj(n-^EI^~V,{N-^E) = ttA^-IV.^Px, (U) 

where tt is defined as tt = (/> — N^Vjcj)- We decompose the shift vector A^' into scalar and intrinsic vector parts as 
Ni = Ni + diip, where diN^ = 0, here the index is raised with 6ij. 

On the comoving time-slices, the scalar field fluctuations vanish, 6(j) — 0, and the three-dimensional spatial metric 
hij is perturbed as (ssj 

=a^e^^6,j, (12) 

where C denotes the curvature perturbation on comoving slices and tensor perturbations have been neglected because 
they do not contribute to the tree-level scalar bispectrum. 
We expand N and A^' in powers of the perturbation ( as 

N = l + ai+a2 + --- , TV, = iv/^^ 7V/^' • • • , ^' = V--! + i/'2 + • • • , (13) 

~ (") 

where «„, Ni and tpn are of order C". 

Now, the strategy is to solve the constraint equations for the lapse function and shift vector in terms of ( and then 
plug in the solutions in the expanded action up to the desired order. It turns out that even for the third order action 
one does not need to use the explicit solution for the constraints past first order 0, ■ 

At first order in C, the solution for equations (jlip with a particular choice of boundary conditions at spatial infinity 

is I, m m 



ai = |:, iv/'^=0, ^, = -^+x, d\^a'^C. (14) 



The second order action is then 



S2^ J dtd^x^ (C^ - a-^c^idCf) . (15) 

At leading order in slow-roll, H is constant and the scale factor may be approximated by that of a pure de Sitter 
universe a = —1/{Ht) and —Ht = ln(— i/r), where r denotes conformal time. The solution of the equation of motion 
derived from the previous action at leading order in the slow- variation parameters and in Fourier space is Q 

C(r,k) = -^^=(l + ^fcc.r)e-'='==^ (16) 

To quantize the curvature perturbation, we follow the standard procedure in quantum field theory. We promote C 
to an operator that is expanded in terms of creation and annihilation operators and mode functions as 

C(t, k) = C(t, k)a(k) + C*(t, -k)at(_k), (17) 

where a and a' satisfy the usual commutation relation [a(k), a^(k')] = {2tt)^ 6^'^'^ (k — k'). 
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IV. THE THIRD ORDER ACTION AND THE BISPECTRUM 



In this section, we will present the standard equations needed in the calculation of the bispectruni of the primordial 
curvature perturbation using the so-called in-in formalism [65l.l66j. In subsection llV Al we shall consider the calculation 
using the well-known field redefinition prescription and in subsection II V Bl we will show that we can obtain the same 
result if we include the contribution from the boundary terms that appear when we simplify the action to the form 
(HH]) by using integrations by parts. 

In order to use the machinery of the in-in formalism to compute the tree-level three-point correlation function (or 
bispectrum) one needs to calculate the cubic-order interaction Hamiltonian, see for example [11] for a review about 
this procedure. 



The field redefinition 



The third order action, ignoring the many boundary terms, has been known since the seminal work by Maldacena [5S 
(for the standard kinetic term case) and can be also found in 0,163 (for the model ©), it reads 



Reduced 



dt (TX 



«Ms(l-^)+2A)|, + ^(.-3-^3c^)CC^ + 5(.-2,s-Hl-c^)C(aC)^ 



a?e d f rj 



2a-C{dC){dx) + ( 3 ) C\ + -{dO{dx)d'x + -{d'Oidx? + 2/(C) 



where we should note that no slow-roll approximation has been made and we define 



,(18) 



(19) 



1 



4a^i/' 



[-{dQ){dO + d-\d,d,{d.,cd,0)\ 



2a?H 



[{dc){dx) - d-\d,d,md,x))] 



(20) 



In Eqs. dUl) and denotes the inverse Laplacian and 5L/S(^\i denotes the variation of the quadratic action 

with respect to the perturbation (. The last term in Eq. (fT8| can be absorbed by a field redefinition of ^ [3], 



C ^ Cn + /(CO- 



(21) 



After this redefinition of variables, there appears a contribution to the interaction Hamiltonian at cubic order coming 
from the originally second-order action (See Eq. ((3T|) below). This term is identical to the symmetric of the last term 
in Eq. ()18p up to a surface term arising from the integration by parts. Neglecting such a surface term, the interaction 
Hamiltonian in conformal time becomes 



H,nt{T) 



E(l-l)+2A)g + 5(e-3 + 3c?)C„C + f (.^2,s + l-c?)C„(5Cn)^ 



-2^C{dQn){dXn) + ^(\ 

ci 2c: \ci 



CCn + f (5Cn)(aX„)a2x„ + ^{d'CnKdXnf 

2a 4a 



(22) 



where prime denotes derivative with respect to conformal time and Xn — o-^Cg ^C4- 
The field redefinition (1211) introduces extra terms in the three-point function as 



(C(X1)C(X2)C(X3)) = (C„(X1)C„(X2)C„(X3)) 



+ ^(C„(xi)C„(x2))(C„(xi)C„(x3)) -Ksym-|-0(ry2c;4(p^)3)^ 



(23) 



where "sym" denotes two terms that result from the preceding one by symmetrizing with respect to Xi, X2 and X3. 
The field redefinition (1211) includes several other terms, however in the previous expression we only displayed the 
contribution of the first term in ((20)) . This is because the omitted terms involve at least one derivative of C and they 
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should vanish when evahiated outside the horizon giving a neghgible contribution to (j23|) . For the reason mentioned 
above, in this paper we approximate the function /(C) as 

m - i^f- (24) 

In the expression p3p , the slow roll parameter 77 has to be evaluated at the end of inflation. This means that 77 might 
become large depending on how inflation ends which in turn would imply that the expansion in terms of 77 would 
cease to make sense. However in most cases one can safely ignore the last term in Eq. (j23p . which is higher order in 
77 because this is also greatly suppressed by powers of the spectrum which is measured to be of order 10^ 
The tree-level three-point correlation function at the time after horizon exit is 

(r!|Cn(Te,ki)C„(re,k2)C«(re,k3)|r^) = / " d™(0| [(„ (Xe , ki )(„ (t^ , ka )C« (Te , kg) , (f )] |0) , (25) 

where jfi) and |0) denotes the interacting vacuum and the free theory vacuum respectively. [, ] denotes the standard 
commutator and the interaction Hamiltonian Hint is used to evolve the free theory vacuum to the interaction vacuum 
at the time when the three-point function is evaluated [58| . 

In Maldacena's calculation [58j and in the following ones [3, H^l the last three terms in Eq. ([2^ which are higher- 
order in the slow-roll expansion were properly neglected because these authors work at leading order in slow-roll. In 
this work, we are only interested in the bispectrum produced by the field redefinition, i.e. the second line of Eq. (|23p . 
In Fourier space it reads [3] 

imoMKiOMKiOMm = (27r)3^(^)(K) ^g^^ %k!h)? ' ^2^) 

where K is defined as K = ki + k2 + kg. 



B. The boundary terms 



In the previous expression for the third order action ()18p we omitted both time and spatial boundary terms as in 
previous calculations. However the total action includes them and in this subsection we will obtain these interactions 
explicitly. In the following, total spatial derivative terms will be omitted because they do not contribute to the three 
point function. It is easy to see that the interaction Hamiltonian of these terms is proportional to K = ki + k2 + kg 
which has to vanish because of momentum conservation imposed by the overall Dirac delta function. 

The third order action without neglecting the boundary terms is given by 

gTotal ^Reduced _j_ ^Boundary (27) 

and the explicit form of the boundary term is 



^Boundary ^ I ^3 ^ ± 
^ ' dt 



'-9a^He + ^adc?-^idcrd\ 

fn^ ■ 1 



(28) 



2c1 " 2aH 

where once again total spatial derivative terms were omitted. One can calculate the three point functions coming from 
the individual terms in this expression, directly using C. The difference in the third order action compared with the 
case of Cn is the last term in (IT5|) and g^°'^"-'^°-^y ^ xhe last term in (US]), 2f{()(6L/dC)\i does not give any contribution 
because the factor {5L/S()\i vanishes after substitution of the mode function. The contributions from the different 
terms in 5*3 °"" "''^^ all vanish except for the one coming from the first term in the third line (using the leading order 
mode function and taking the limit Tg — > 0). 

The contribution from the first term in the third line, i.e. 
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can be calculated using the in-in formalism to find 

(f7|C(Te,ki)C(re,k2)C(re,k3)|r!) 



IjecJ^a^iT,, ki)C(Te, k2)C(Te, k3)C*(Te, ki)C(Te, k2)C * {t,, kg 



5 perms. 



re->0 



(30) 



where we used the mode function solution (ITCl) . This bispectrum is of the so-called local type and exactly agrees with 
Eq. (|26l) . This result shows explicitly that the total time derivative terms cannot a priori be neglected [53, [S^ • The to- 
tal bispectrum is the sum of the expression (pO| (or equivalently (|26)) ) with the bispectrum (Cn(0, ki)(^„(0, k2)Cn(0, kg)), 
which is calculated using Eq. ([25]) . This is the result that Chen et al. Q got. They performed the calculation by the 
field redefinition prescription and used Wick's theorem as described in the preceding subsection. 

In the approach of using the field redefinition explained in the preceding section, we did not take into account the 
boundary action g^°^^'^'^^y ^ jf -^g perform the change of variables C = Cn + fiCn) without neglecting the surface term, 
it generates a third order action from the second order action (jl5p as 



dtd'xi^2)fiC,,) ^ 







+ 


/ dtd^x^- ( 


1 


J dt \ 


+ 


f dtd^x— ( 


1 


J dt \ 



€ 

3- 



o 4 SnVl I : 



(31) 



where again total spatial derivative terms were omitted because they do not contribute to the three point function, 
and the approximate equality "sa" means the truncation given in Eq. (1241) . e.g. we neglect total time derivatives terms 
that do not contribute to the leading order three point function. The first and second terms, respectively eliminate 
the last terms in 5'|?ed«ce(i pj-oportional to the linear equations of motion and the boundary term that contributes to 
the bispectrum. 

The boundary terms (j28p are necessary to erase from the action the terms with second order time derivative on ^ 
contained in the last interaction in Eq. ()18|) . Without these boundary terms, the second order time derivative terms 
can remain in the action. These problematic interactions from the point of view of the Hamiltonian formalism are 
not present in the original action Eq. and they are not generated by the procedure of replacing the solution of 
the constraints back into the action. In fact, these terms are generated by the integrations by parts, and therefore 
the inclusion of the boundary terms naturally takes care of them. 

The calculations presented above clearly show that the boundary terms can affect the result of the computa- 
tion. However, if the results of the computation depend on the choice of the total derivative terms in the action 
^ dt d'^ X dB{C,, C,) / dt, one might be worried of how one can correctly choose B. The terms in B containing C, produce 
the second order time derivatives in the action. As mentioned above, such problematic terms should vanish in total. 
This condition completely determines the terms containing C, in B. In contrast, there is no criteria to choose the terms 
that do not contain C,. However, those terms do not contribute to the equal-time expectation values of ^s. In the 
path integral expression for the expectation value of an operator 0{t), adding such boundary terms given by S(C) is 
equivalent to replacing the operator O with O' = e"'-^'' ^B{Qq^'i f d xB{Q ^ However, since B{() does not contain the 
conjugate momentum of C, B commutes with O. Therefore C reduces to O and hence the results are not affected by 
such boundary terms. 



V. CONCLUSION 



We have computed the total time derivative interactions in the third order action for the comoving curvature 
perturbation in general k-inflation. In previous calculations of the bispectrum these boundary terms have been 
ignored. In this note we have shown explicitly that a priori they are important and should not be neglected freely. 
These boundary interactions are necessary to erase the terms with second-order time derivatives on C in the action, 
generated by the integrations by parts. Total spatial derivative terms can be safely ignored because their contribution 
for the bispectrum is proportional to the sum of the three momentum vectors which has to be zero due to momentum 
conservation. 

From all the boundary terms that appear after many integrations by parts in the action only one of them gives a 
non-zero bispectrum at leading order in the slow-variation expansion. We have shown that the bispectrum produced 
by this term is equal to the bispectrum produced in the usual field redefinition prescription, thus our results agree 
with previous results in the literature. 
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The main conclusion of this note is that in the calculation of the bispectrum in general k-inflation one can ignore 
all the boundary terms that appear when one simplifies the action but then one has to perform a field redefinition 
to eliminate terms in the action that are proportional to the first order equation of motion. On the other hand, one 
might choose to keep all the boundary terms and calculate the bispectrum using the usual method without the need 
to do the field redefinition. We have shown that in the end the bispectrum of the curvature perturbation is the same 
in both procedures. 

An important lesson in computing the tree-level bispectrum is that we should basically use the reduced action 
written in terms of the physical variables that does not contain second order time derivatives in total. If second order 
time derivatives are contained, they must be eliminated by integration by parts. Only when the action takes the 
canonical form without second order time derivatives, one can use the expression for the third order action as it is in 
the path integral expression. 

Note added: While we were writing up this work, the paper (67| appeared in the arXiv. They cite a paper 
in preparation [68j where the authors also argue that the inclusion of the boundary terms accounts for the terms 
introduced by the field redefinition. 

Note added in version 2: Soon after this paper appeared on the arXiv, Ref. [68*1 appeared as [6^. Their sec- 
tion 3.1.2 contains similar, but independent, arguments on the role of the boundary terms. Some time after that, 
Ref. [70| appeared on-line and it also discusses the relation between field redefinitions, boundary terms and gauge 
transformations in the computation of the bispectrum. 
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